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Abstract 

The influence of quantum fluctuations on electron transport through small 
metallic islands with Coulomb blockade effects is studied beyond the pertur- 
bative regime. In tunnel junctions with low resistance higher order coherent 
processes and "inelastic resonant tunneling" become important. We present 
a path integral real-time description, which allows a systematic diagrammatic 
classification of these processes. Quantum fluctuations renormalize system 
parameters and lead to finite lifetime broadening. Both effects are observ- 
able in the gate voltage dependence of the nonlinear conductance. The finite 
bias voltage introduces an energy scale up to which quantum fluctuations are 
probed. It can be larger than the temperature, which makes quantum effects 

experimentally observable above thermal fluctuations. 
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Electron transport through mesoscopic metallic islands is strongly influenced by the 
charging energy associated with the low capacitance of the junctions A variety of 

single-electron phenomena, including Coulomb blockade at low temperatures and Coulomb 
oscillations of the conductance as a function of a gate voltage have been observed. A master 
equation description of sequential tunneling |T[ is sufficient as long as the conductance of 
the barriers is low, «o = h/(4ir 2 e 2 RT) <C 1. In general, quantum fluctuations and higher 
order coherent tunneling processes should be considered |§-f|. This includes "cotunneling" 
ill , where in a second order coherent process electrons tunnel via a virtual state of the 



island. Furthermore, resonant tunneling plays a role. In comparison to the well known 
case of independent electrons we encounter here two complications: (i) The metallic system 
contains many electrons, and with overwhelming probability different electron states are 
involved in the different transitions of the coherent process (this is denoted by "inelastic"), 
(ii) The strong Coulomb interaction cannot be accounted for perturbatively. 

In the present article we develop a systematic diagrammatic technique to identify the 
processes of sequential tunneling, cotunneling and inelastic resonant tunneling. We study 
the time evolution of the density matrix. In an earlier paper || we had formulated the 
problem, after a separation of charge and fermionic degrees of freedom, in a many-body 
expansion technique. Here we reformulate it in a real-time path-integral representation JIT] . 



The latter is familiar from the work of Caldeira and Leggett ]12| and Feynman and Vernon 



13fl who studied dissipation in quantum mechanics. Dissipation associated with tunneling 



of electrons was investigated in Refs. |fl^j3[1. Similar problems arise in the context of local, 



strongly correlated Fermi systems like the Kondo and Anderson model [|15[]. An essential 
step in the present work is a transformation of the path-integral description of electron 
tunneling from a phase to a charge representation 0,0. 

As examples we study the single electron transistor || and the electron box. In the 
transistor a metallic island is coupled by two tunnel junctions (L, R) to two electrodes, 
and further capacitively to a gate voltage Vq. A transport voltage V = Vi — Vr drives a 
current. The charging energy is E c h(n) = (ne — Qg) 2 /{2C). It depends on the number 



of excess electrons n on the island and on the continuously varying external charge Q G = 
C G V G + C L V L + C R V R . The total island capacitance C = C L + C R + C G defines the scale of 
the charging energy, E G = e 2 /2C. The electron box, consisting of an island coupled via a 
tunnel junction to one electrode and capacitively to the gate voltage, is described similarly. 

At low temperature the electron number n{Q G ) in the electron box increases in unit 
steps by tunneling as Q G is increased. In the transistor tunneling occurs in lowest order 
perturbation theory only if the electrochemical potential of one electrode is high enough to 
allow an electron to enter the island, say eV L > E ch (n+1) — E ch (n), while the electrochemical 
potential of the other electrode allows the tunneling process to that electrode, i.e. E c h(n + 
1) — E ch (n) > eV R . Within the window set by these conditions the current is AR t I(V, Q G ) = 
V - 4[Q G - (n + l/2)e] 2 /(C 2 V). I.e. the nonlinear conductance, G(V, Q G ) = dl/dV, shows 
as a function of Q G an e-periodic series of structures of width CV with, at T = 0, vertical 
steps at its edges. At finite temperature (n(Q G )) and the steps are washed out. 

Quantum fluctuations further wash out these steps. Below we describe these processes 
diagrammatically and determine their effect on the nonlinear conductance. In the most 
interesting case we can re-sum the diagrams and obtain closed expressions for the stationary 
density matrix and the spectral density describing the charge excitations of the system. 
Our main findings are: (i) Quantum fluctuations renormalize the energy and conductance, 
(ii) The coherent resonant tunneling processes give rise to a broadening of the charge state 
levels. The effects are observable experimentally in the nonlinear conductance. Here the 
finite bias voltage introduces an energy scale, up to which quantum fluctuations are probed. 
Since eV can easily be chosen larger than the temperature, quantum effects are not hidden 
below thermal broadening and become observable in a realistic experiment. 

The description of the single electron transistor is based on the Hamiltonian H = + 
H R + Hj + H ch + H t L + H t L . The first terms describe noninteracting electrons in the metallic 
left and right lead and island. These are treated as reservoirs, i.e. their electron distributions 
are assumed to be thermal with electrochemical potential which depends on the state. The 
Coulomb interaction H c h depends on the charge on the island, as expressed by the charging 
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energy E c h(n) given above. Charge transfer is described by standard tunneling Hamiltonians 
Ht, r = 12k&r,q£i ,aT r c\ a c qa + h.c. , where r = L,R. The matrix elements and densities of states 
determine the tunnel conductances, R^ 1 = (ine 2 /h)N r (0)Nj(0)\T r \ 2 . We consider "wide" 
metallic junctions with many transverse channels N c h ^> 1. Hence "inelastic" higher order 
tunneling processes, involving different electron states for each step, dominate over "elastic" 
processes which involve the same state repeatedly. 

We describe the system of interacting electrons in a path integral representation. A 
Hubbard-Stratonovich transformation, used to handle the capacitive interaction, introduces 
as collective variable the phase (p, the quantum mechanical conjugates of the charge on 
the island. The phases in the electrodes are fixed by the applied voltages, (p r = eV r t. 
Next, the electronic degrees of freedom can be traced out, followed by an expansion of the 
electron propagators. Since we consider wide junctions, only simple loops need to be retained 
||16|| . Their iteration introduces in each order a factor N c h, hence they dominate over more 
complicated higher order loops. After this stage the system is described by a reduced density 
matrix p{{fi},{f2})- It can be expressed by an effective action, which depends on the 
phases <p a corresponding to the forward and backward propagator o = 1,2. The structure 
of the theory is familiar from Refs. P~5|,[T2[] , where a quantum system coupled to a harmonic 



oscillator bath has been considered, and from Refs. ]14]|§, where electron tunneling has been 
described. The effective action contains the charging energy of the system appropriate to 
E c h. The tunneling couples the forward and backward time propagators. For each junction 
we have [0,0] 

S t ,r[<PiM=±™ ]T [ tf dt f dt'a a /{t-t') cosf^Ct) - <p a ,(t')} ■ (1) 

a,a> = 1,2 l ' li 

The kernels are given in Fourier space by a^iu) = (— l) a+1 a~(u>), ot° ,2 {uj) = (— l) ff a^(u) 
with af{u) = ±a 0tr (lu — eV r ) [exp(±/3(u; — eV r ) — l] -1 . 

An important step for a systematic description of tunneling processes is the change from 
the phase to a charge representation ||, accomplished by (h — k — 1) 

p(t f ;n lf ,n 2f ) = 



J2 p(ti;n u ,n 2 i) [ d(p lf dip 2 fdip li dip2i f ' Vip^U) f ^ T><p e (\J) f X>\«,(U) f V\ e (U) 
exp ji^£ t)r [<pi,<p 2 ] + Yl (-l) a (^ n <ri¥cri - inafCfcrf + iSch[n a ) - i J dtn^^^j . 



(2) 



The charging energy is S c h\n c \ = f t f dt^{n a e — Qg) 2 ■ In systems with discrete charges the 
integrations include a summation over winding numbers 0. Next we expand the tunneling 
terms exp(iSt) and integrate over ip a . Each of the exponentials exp[±z<^o-(£)] describes 
tunneling of an electron at time t on the forward or backward branch, a = 1 or 2, respectively. 
These changes occur in pairs in each junction r = L, R and are connected by a%' a (t — t'). 
Each term of the expansion can be visualized by a diagram. Several examples are displayed 
in Fig. ([!]). There is a closed time-path consisting of two horizontal lines, corresponding to 
the forward and backward propagator between £j to tf. Along the time-path vertices are 
arranged, connected in pairs by (dashed) tunneling lines, either within one propagator or 
between the two propagators. 

We start from a density matrix which is diagonal, p(ti) = p(°\n)\n)(n\. During tran- 
sitions the system is in an off-diagonal state. We denote the sum of all diagrams, starting 
and ending in the diagonal states n and n', respectively, by H n ,n'- It can be expressed by an 
irreducible self-energy S n n /, which describes the transitions, in the style of a Dyson equation 

n ni n' = H^n + ^2 n-n,n" £n",n' ^°'n> . (3) 
n" 

Here IT^ is the free propagator. We identify the solution of Eq. (|3|), multiplied with 
P(°)(n), as the stationary distribution function Xm P i n )^n,n' — P st (n'). The sum rule 
J2 n ' S n , n ' = (see Ref. |TT| for further details) implies 



= -P s \n) ]T £ ft , n , + J2 P St (n')Z n >,n ■ (4) 

I.e., we recover a stationary master equation. In general the irreducible self-energy S n / jn 
yields the transition rates of all possible correlated tunneling processes. We, furthermore, 
note that the stationary solution P st {n) is independent of the initial distribution P^(n). 
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For illustration we evaluate all diagrams which contain no overlapping tunneling lines, 
as visualized on the left hand side of Fig. ([!]). In this case the irreducible self-energy parts 

E« ± i = 2ma ± (±AE±) , E« = -2ni £ a ± (±AE^ h ) , (5) 

± 

where AEf h = E^nil) — E c h(n), reproduce the well-known single electron tunneling rates. 
In situations where single electron tunneling is suppressed by Coulomb blockade the lowest 
order contribution to the current arises due to cotunneling. It is described by a diagram, 
also shown in Fig. ([!]), with tunneling processes in the left and in the right junction, where 



the corresponding lines — t' L ) and or^h — t' R ) overlap in time. Applying our rules JTT 

we reproduce the well-known cotunneling rate. 

The perturbative approach is sufficient for a o m (^0 < 1. At larger values of the 
conductance resonant tunneling processes get important. To proceed we have to find a 
systematic criterion which diagrams to retain. For this we note that during a tunneling 
process the reservoirs contain an electron excitation. Our criterion is to take into account 
only matrix elements of the density matrix which differ at most by two excitations in the 
leads or (equivalently) in the island. This means that in the diagrams any vertical line cuts 
at most two tunneling lines. We, furthermore, concentrate on situations where only two 
charge states, n = 0,1, need to be considered. This is sufficient when the energy difference 
of the two states A = E ch {l) — E ch {0) and the bias voltage eV = eV^—eV^ are low compared 
to the energy, Ec, associated with transitions to higher states. The combination of both 
restrictions implies that the diagrams contain no crossing tunneling lines, which allows us 
to evaluate the irreducible self-energy analytically. 

Using the notations a r (u) = a^{uj) + a~(ui), ai(u) = J2r oz r (uj) and ao = J2 r a o,r we find 

E 0) i = -£ ,o = 2Tci ~x~ > Sl >° = "^i- 1 = 27 ™~\~ ^ 
with A± = y'dcja ± (cj)|7r(cj)| 2 , X = J did \tt(lu)\ 2 , (7) 

and tt(cj) = — ^ , a(u) = - [ dJ — ^— ^ . (8) 

Inserting these quantities in Eq. @ we arrive at the stationary, normalized probabilities 
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P,f = A_ and Pf = A+ , with A+ + A_ = 1. 

The expression for the current at time t in the junction r can be written as 

I r (t) = 4vue /* dt'J2^ r ' a (t - t')(smbi(t) - MOD , (9) 

J-oo a 

where the expectation value is taken with the density matrix discussed above, with t = tf. 
We, therefore, study the correlation functions describing charge transfer at times t and t' 

C>(t,t') = -z( e -^(V^')) , C<(t,t') = i^We-**®) . (10) 

We further introduce a spectral density for charge excitations A(v) = j^[C < (ou) — C > (u;)]. 
Within the approximations described above we find 

C>H = ±2vr^a r H/[±(cu-eK)]kM| 2 , (11) 

r 

A{u) = a(uj)\TT{uj)\ 2 (12) 

and I r = W I du>Y, ^ U) ?} U) A(u)\J{u - eV r .) - f(u - eV r )\ . (13) 

These results satisfy conservation laws and sum rules, J2r Ir = and / duA(u) = 1, and the 
equilibrium relations between correlation functions and the spectral density. The classical 
result is recovered in lowest order in a , where A^(lj) = 5(u — A ). Quantum fluctuations 
yield energy renormalization and broadening effects, which enter in the spectral density via 
the complex self-energy a{uS) given in Eq. @. 

In equilibrium, for Vr = Vl ) the transistor is equivalent to the single electron box. The 
average electron number becomes (n) = J dujf(uj)A(uj). In the classical limit it reduces to 
(n d ) = /(A ), where the energy difference A = E c (l — ICgVgI^ depends on the gate 
voltage. It shows a step at Qg = e/2, which is smeared by temperature. At larger values 
of a or lower temperature we have to include the self-energy <r(u) in the spectral density. 
The limits T, V = and |u;| < T,V can be analyzed analytically. In the former, 

A{u) - M . AM5M , (14) 

V ' A [u - A{uo)f + [ttA(cu)«(cu)] 2 V ; 

= l + 2a ln(§) ' l + vr^M 2 ' (15) 

ol{uj) = — --b— . (16) 

1 ' l + 2a ln(§) 1 ' 



The spectral density A(u) has a maximum at the renormalized energy difference A, obtained 
from the self-consistent solution of A = A (A) [§,0. It further is broadened by nAa, where 
a = a(A). Due to quantum fluctuations the step of the average charge (n) in the electron box 
at the degeneracy point Ao = is washed out even for T = 0. At finite T the slope at Ao = 
shows an anomalous temperature dependence, d(n) /<9A |a =o — — {4T[l + 2a ln (j-y)] 2 } -1 - 
A pronounced and experimentally accessible signature of quantum fluctuations is con- 
tained in the nonlinear response of a transistor. We study the differential conductance 
G(V, Qg) = dI Qy S> at finite voltages. In this case the spectral density in the limit 
\u\ < T, V is involved. Here we find a renormalization of A and «o by a factor Z with 
Z^ 1 = 1 + 2aoln(Ec/ max{el^/2, 2ttT}) and a broadening of the spectral density which is 
given by ZaiuS). Now eV provides an energy scale, and renormalization and life-time effects 



are probed over a wide energy range even at zero temperature. The result of Eq. (|13|) is 
plotted in Fig.(fJ). For comparison we also show the classical result, where the conductance 
is nonzero only in the range |A | < with vertical steps at the edges. Fig.(|2|) clearly dis- 
plays the renormalization effects and, moreover, the finite life-time broadening. The width 
of the structure as a function of Aq/V is Z~ x which depends logarithmically on V. These 
effects are observable in an experiment with realistic parameters. In particular, a transport 
voltage can be applied which exceeds the temperature, kT <C eV < Eq- Hence thermal 
smearing is still weak, and quantum effects are visible. 

In conclusion, we have presented a systematic description of coherent single electron 
tunneling processes including resonant tunneling in a transistor. We propose to measure 
the nonlinear conductance in order to study renormalization effects and lifetime broadening 
induced by these quantum effects. 

We acknowledge the discussions with G. Falci, A.D. Zaikin and G. Zimanyi. This work 
was supported by the 'Sonderforschungsbereich 195' of the DFG and the Swiss National 
Science Foundation (H.S.). 
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FIGURES 

FIG. 1. Example of a diagram showing various tunneling processes: on the left sequential 
tunneling in the left and right junctions, then a term which preserves the norm, next a cotunneling 
process, and on the right resonant tunneling. 

FIG. 2. The differential conductance at T = in the nonlinear response regime as function of 
the gap energy normalized to the transport voltage V. We consider a symmetric bias and choose 
«o,L = a ,R = 0.05, and (1) eV/E c = 0.1 , (2) eV/E c = 0.01 , (3) eV/E c = 0.001 . For comparison, 
(0) shows the result obtained in lowest order perturbation theory. 
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